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Abstract 

In this paper we consider the design of spectrally efficient time-limited pulses for ultra-wideband (UWB) systems 
using an overlapping pulse position modulation scheme. For this we investigate an orthogonalization method, which 
was developed in 1950 by Lowdin l!l,|2[]. Our objective is to obtain a set of N orthogonal (Lowdin) pulses, which 
remain time-limited and spectrally efficient for UWB systems, from a set of N equidistant translates of a time-limited 
optimal spectral designed UWB pulse. We derive an approximate Lowdin orthogonalization (ALO) by using circulant 
approximations for the Gram matrix to obtain a practical filter implementation as a tapped-delay-line [7]- We show 
that the centered ALO and Lowdin pulses converge pointwise to the same square-root Nyquist pulse as N tends to 
infinity. The set of translates of the square-root Nyquist pulse forms an orthonormal basis for the shift-invariant-space 
generated by the initial spectral optimal pulse. The ALO transformation provides a closed-form approximation of 
the Lowdin transformation, which can be implemented in an analog fashion without the need of analog to digital 
conversions. Furthermore, we investigate the interplay between the optimization and the orthogonalization procedure 
by using methods from the theory of shift-invariant-spaces. Finally we relate our results to wavelet and frame theory. 
Keywords: Lowdin Transformation, Shift-invariant spaces, PPM UWB transmission, FIR filter design, canonical 
tight frame 



1. Introduction 

We consider in this work high data rate transmission in the ultra wideband (UWB) regime. To prevent disturbance 
of existing systems, e.g. GPS and UMTS, the Federal Communications Commission (FCC) released [5] a very low 
power spectral density (PSD) mask for ultra-wideband (UWB) systems. To ensure that sufficiently high signal-to- 
noise ratio (SNR) is maintained in the frequency band F = [0, 14]GHz, as required by the FCC, the pulses have 
to be designed for a high efficient frequency utilisation. This utilisation can be expressed by the pulses normalized 

In- 



effective signal power (NESP) 



Several pulse shaping methods for pulse amplitude and pulse position modulation 



(PAM and PPM) were developed in the last decade based on a FIR prefiltering of a fixed basic pulse. Due to the high 
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sampling rates in UWB this FIR filtering is realized by a tapped delay line |9J]. A SNR optimization under the FCC 
mask constraints then reduces to a FIR filter optimization 16j-|9J]. 

Since the SNR is limited, the amount of signals can be increased to achieve higher data rates or to enable multi-user 
capabilities. For coherent and synchronized transmission over memoryless AWGN channels, an increased number N 
of mutually orthogonal UWB signals inside the same time slot, known as A^-ary orthogonal signal design, improves 



the BER performance over Eb/No and hence the achievable rate of the system [ 10, Ch.4]. 



Combining spectral shaping and orthogonalization is an inherently difficult problem being neither linear nor con- 
vex. Therefore most methods approach this problem sequentially, e.g. combining spectral optimization with a 
Gram- Schmidt construction [9, nj I2L This usually results in an unacceptable loss in the NESP value of the 
pulses J13I [14I . Moreover these orthogonal pulses are different in shape and therefore not useful for PPM. Further- 
more, a big challenge in UWB impulse radio (UWB-IR) implementation are high rate sampling operations. Therefore, 



an analog transmission scheme is desirable [ 14- 16] to avoid high sampling rates in AD/DA conversion 11511 . 



Usually, PPM is referred to an orthogonal (non-overlapping) pulse modulation scheme. To achieve higher data 
rates in PPM, pulse overlapping was already investigated in optical communication I17I and called OPPM. An ap- 
plication to UWB was studied for the binary case with a Gaussian monocycle [18]. To the authors' knowledge, 
no orthogonal overlapping PPM (OOPPM) signaling has been considered based on strictly time-limited pulses. In 
this work we propose a new analog pulse shape design for UWB-IR to enable an almost OOPPM signaling which 
approximate the OOPPM scheme up to a desired accuracy. 

In our approach, we first design a time-limited spectral optimized pulse p and perform afterwards a Ltiwdin 
orthogonalization of the set of 2M+1 integer translates \p(- -k)}f = M , which span the function space V M (p). This 
orthogonalization method provides an implementable and stable approximation p° M to a normalized square-wot 
Nyquist pulse p° (all integer translates are mutually orthonormal) . Using the Fourier transformation for p given by 
p(v) = L p(t)e~' 2 "' v dt for every v e R, the square-root Nyquist pulse p° can be expressed in the frequency-domain 

P (v) = - — (1) 

^ k \P(y-k)\ 2 

for v almost everywhere, which is known as the orthogonalization trick JjjJ]. 

Usually in digital signal processing, see Fig.Q] the time-continuous signals (pulses) in L 2 are sampled by an analog 
to digital (AD) operation (either in time or frequency) to obtain time-discrete signals in C 2M+1 . Then a discretization 
of the orthogonalization in ([TJ yields a finite digital transformation D M to construct a discrete signal in the frequency 
domain which has again to be transformed by a DA operation to obtain finally an approximation of the square-root 



Nyquist pulse |2C ] 



Instead of using such an AD/DA conversion to operate in a discrete domain (depictured in Fig.[TJby the digital box), 
we use the democratic Lowdin orthogonalization^ B M (in the analog box), found by Per-Olov Lowdin in | lj 31, 



2 In fact, the Lowdin orthogonalization and the orthogonalization trick are both orthonormalizations, but for historical reasons we will refer them 
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Figure 1 : Analog and discrete approximation methods in time and frequency domains 



where all 2M + 1 linear independent pulse translates are involved simultaneously by a linear combination to generate 
a set of time-limited mutual orthonormal pulses. Hence the Lowdin orthogonalization is order independent. The 
Ltiwdin pulses constitute then an orthonormal basis for the span V M (p) of the initial basis \p{- - k)}f = M . Moreover, 
as we will show in our main Theorem [2] the Lowdin orthogonalization B M is a stable approximation method to the 
construction in ([TJ and operates completely in the analog domain. Note, that a discretization of ([]]) generates neither 
shift-orthonormal pulses nor a set of mutual orthonormal pulses. The important property of the Lowdin method is its 
minimal summed energy distortion to the initial basis. It turns out that all orthonormal pulses maintain the spectral 
efficiency "quite well". 

As M tends to infinity the Lowdin orthogonalization B = B°° applied to the initial pulse p delivers the square-root 
Nyquist pulse p°, which allows a real-time OOPPM system with only one single matched filter at the receiver. Since 
the Lowdin transform B M is hard to compute and to control we introduce an approximate Lowdin orthogonalization 
(ALO) B M and investigate its stability and convergence properties. It turns out that for fixed M the transformations 
B M and B M are both implementable by a FIR filter bank (realised as a tapped delay line) like the spectral optimization 
in |9|. We call B M p and B M p approximate square-root Nyquist pulses, since we observed that even for finite M our 
analog approximation yields time-limited pulses with almost shift-orthogonal character since the sample values of 
the autocorrelation are below a measurable magnitude for the correlator. Hence such a construction of approximate 
square-root Nyquist pulses seems to be promising for an OOPPM system. 

The structure of this paper is as follows: In Section [2] we introduce the signal model and motivate our spectrally 
efficient A^-ary orthogonal overlapping PPM design for UWB systems. Section[3]presents the state of the art in FCC 
optimal pulse shaping for UWB-IR based on PPM or PAM transmission by FIR prefiltering of Gaussian monocycles 
which is a necessary prerequisite for our design. To develop our approximation and convergence results in Section |4] 
we introduce the theory of shift-invariant spaces in Section |4~T1 and in Section l4~2l the Lowdin orthogonalization for 
a set of N translates to provide a N-axy orthogonal overlapping transmission. Our main result is given in Section l43l 
where we consider the stability of the Lowdin orthogonalization B M (for M increasing to infinity) and develop for 



as orthogonalization methods. 
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this a simplified approximation method B M , called ALO. In Section [5] we study certain properties of our filter design 
and investigate the combination of both approaches. Furthermore, in Section [57X1 we develop a connection between 
our result and the canonical tight frame construction. Finally in Section[6] we demonstrate that the ALO and Lowdin 
transforms yields for sufficiently large filter orders compactly supported approximate square-root Nyquist pulses, 
which can be used for OOPPM having high spectral efficiency in the FCC region. Moreover, the Lowdin pulses 
provides also a spectrally efficient (2M + l)-ary orthogonal pulse shape modulation (PSM) Oil . 



2. Signal Model 

In this work we will consider finite energy pulses p, i.e. the set L 2 :- { p : 

. n 1 

integrable functions with norm ||-|| L 2 := V(-, ■) induced by the inner product] 

(p,q):= I p(t)q(i)dt. 

Jr 

To control signal power in time or frequency locally we need pulses in L°° :- J p 
essentially bounded, i.e. functions with a finite L°°-norm, given as 



j| L 2 < oo! of square- 



(2) 



-> J which are 



\\P\L •= esssup{|p(f)|), 

reR 



(3) 



where the essential supremum is defined as the smallest upper bound for \p{t)\ except on a set of measure zero. If the 
pulse is continuous than this implies boundedness everywhere. UWB-IR technology uses ultra short pulses, i.e. strictly 
time-limited pulses with support contained in a finite interval Xci We call such L 2 - functions compactly supported 
in X and denote its closed span by the subspace L 2 (X). The coding of an information sequence \d„} = {d„}„ £ z is 
realized by pulse modulation techniques Jl?! 22I of a fixed normalized basic pulse p e L 2 ([0, T p ]) with duration T p . 

A relevant issue in the UWB-IR framework and in our paper is the spectral shape of the pulse. In this section we 
will therefore summarize the derivation of spectral densities for common UWB modulation schemes such as PAM 



(ZD, PPM ft 



23 



2611 and combinations of both [27] to justify our spectral shaping in the next section. Antipodal 



PAM and iV-ary PPM are linear modulation schemes which map each data symbol d„ to a pulse (symbol) 54(f) with 
the same power spectrum £ |/5(v)| 2 . If we fix the energy £> of the transmitted symbols and the pulse repetition time 
(symbol duration) T s , we will show now for certain discrete random processes ( e.g. i.i.d. processes yfl) that the 
power spectrum density (PSD) of the transmitted signals is given by 



S u (v) 



&\p(v)\ 2 



(4) 



Hence an optimization of the pulse power spectrum to the FCC mask Sf(v) over the band F in Section[3]increases the 
transmit power. To be more precise, PPM produces discrete spectral lines, induced by the periodic pulse repetition, 



The bar denotes the closure for sets and the complex conjugation for vectors or functions. 
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the use of uniformly distributed pseudo-random time hopping (TH) codes c„ e [0, Ay was suggested to reduce this 



effect and to enable multi-user capabilities [23 



25 



28 



29]: 



u(t) = 2^ yfSp(t - nTf - c„T c - di„ /NfS T). 



(5) 



In 12811 this is called framed TH by a random sequence, since the coding is repeated in each frame Nf times with 
a clock rate of 1/7/. Hence NjTj - T s is the symbol duration for transmitting one out of N symbol waveforms 
representing the encoded information symbol d n e {0,...,N - 1). To prevent ISI and collision with other users, the 
maximal PPM shift T and TH shift T c have to fulfill NT < T c and N C T C < Tf. To ensure mutual orthogonality 
of all symbols one requires T > T„ . The PSD for inde pendent discrete i.i.d. processes {c n }, {d„} follows from the 

i QQ tot i 



Wiener-Khintchine Theorem 



S u (v) = 



(5)], 127 



31]. 



(6) 



with |Gjg(v)| 



1 I sin(7TvT c N c ) 



N C N I sm(nvT c ) 



sin(7TvTN) 



and 6 is the Dirac distribution. 



(7) 



sm(nvT) 

However, a more effective and simple reduction method without the use of frame repetition (Nf = 1) or random 
TH (N c = 1) has been proposed in flQ]- Here antipodal PAM with a„ e {-1, 1) is combined with N-stry PPM 
modulation, for NT < T s 



u(t) — Vfi a nP(t - nT s - d„T). 

n 

The PSD for such i.i.d. processes is well known Esec.4.3]: 



(8) 



S«(v) 



£|/3(v)| 2 



l2 |E[a] •E[e-'' 2lrvc,7 'l| _ / 
E y\ - |E [a] ■ E [e-^f + 1 ^ 5 L _ » 



(9) 



since we have Gp(v) = E [a]E |e' 2lrvflrr ] in Q. For an i.i.d. process {a„} with expectation E [a] — and variance 
E |a 2 ] = 1 the PSD reduces to (@}. Hence the effective radiation power is essentially determined by the pulse shape 
times the energy & per symbol duration T s and should be bounded pointwise on F below the FCC mask Sf 



s H (v) = e 



\P(y)\ 2 



< S P (v) for all v e F. 



(10) 



The optimal receiver for A^-ary orthogonal PPM in a memoryless AWGN channel with noise power density No is the 
coherent correlation receiver. The uncoded bit rate Rt, and average symbol error probability P s is given as II Hi 4. 1 .4] 



logA^ 

R b = — — and P S (S) < (N - 1) erfc 




(11) 



where erfc is the complementary error function. Hence, a performance gain for fixed T s is achieved by increasing Af 
and/or &. 
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Increasing N. Usually non-overlapping pulses are necessary in PPM to guarantee orthogonality of the set {p(- - 
riT)} :- \p{- - nT)} ne z of pulse translates, i.e. T > T p . For fixed T p this limits the number of pulses N in [0, T s ] 
and hence the rate In this work we will design an orthogonal overlapping PPM (OOPPM) system by keeping 
all overlapping translates mutually orthogonal. But such square-root Nyquist pulses are in general not time-limited, 
i.e. not compactly supported. In fact, we will show that for a particular class of compactly supported pulses a 
non-overlapping of the translates is necessary to obtain strict shift-orthogonality. However, we derive overlapping 
compactly supported pulses approximating the square-root Nyquist pulse in (Q} and characterize the convergence. 
These approximations to the square-root Nyquist pulse allow a realizable A^-ary OPPM implementation based on FIR 
filtering of time-limited analog pulses. 

Increasing S. The maximization of & with respect to the FCC mask was already studied in 
filtering is used to shape the pulse such that its radiated power spectrum efficiently exploits and strictly respect the 
FCC mask. Note that the FCC regulation in (TlOb is a local constraint and does not force a strict band-limited design, 
however fast frequency decay outside the interval F is desirable for a hardware realization. 



7 ] where a FIR pre- 



Our combined approach now relies on the construction of two prefiltering operations to shape a fixed initial pulse. 
The first filter shapes the pulse to optimally exploit the FCC mask and the second filter generates an approximation to 
the square-root Nyquist pulse. The filter operations can be described as semi-discrete convolutions of pulses p€L 2 
with sequences c e ( 2 := {a | ||a||J 2 =£„ \a„\ 2 <oo) 

p*' T c:=2]c„p(--nT), (12) 

neZ 

at clock rate l/T. If we restrict ourselves to FIR filters of order L the impulse response becomes a sequence c € l 2 (L) 
which can be regarded as a vectors in C L . Hence we refer for c e C L to FIR filters and for c e t 2 to IIR filters. 



3. FCC Optimization of a Single Pulse 

The first prefilter operation generates an optimized FCC pulse p. To generate a time-limited real-valued pulse 

) and 

y, 



we consider a real-valued initial input pulse q e L 2 ([-T q /2, T„/2]) and a real-valued (causal) FIR filter g e R L . A 



see also Section [6] The prefilter operation is 



common UWB pulse is the truncated Gaussian monocycle q ]6, 
then: 

L-\ 

Pit) = (q 4 8)(0 = 2 8kq(t - kT ) (13) 

(t=0 

which results in a maximal duration (support length) T p — (L— 1)Tq + T q of p. 

To maximize the PSD according to (TTOb we have to shape the initial pulse by the filter g to exploit efficiently the 
FCC mask Sp in the passband F p c F, i.e. to maximize the ratio of the pulse power in F p and the maximal power 



allowed by the FCC 

770?):= f \p(v)\ 2 dvl f S F (v)dv. (14) 
Jf p I Jf p 

This is known as the direct maximization of the NESP value J]{p), see J^]. Here we already included the constants & 
and T s in the basic pulse p. If we fix the initial pulse q, the clock rate 1 /To and the filter order L, we get the following 
optimization problem 



max 



v(q*' To &) 



(15) 



s.t. VveF: |g(v)| 2 • \q(v)\ 2 < S F (v), 



where g denotes the I/Tq periodic Fourier series of g, which is defined for an arbitrary sequence c € i 1 as 

oo 

c(v) = c«e~ 2!dmT °- 



(16) 



Since g € ( 2 (L) the sum in (fT6l i becomes finite for g. Liu and Wan 113 311 studied the non-convex optimization problem 
( TT3b with non-linear constraints numerically with fmincon, a Matlab program. The disadvantage of this approach 
lies in the trap of a local optimum, which can only be overcome by an intelligent choice of the start parameters. 

Alternatively dT3T > can be reformulated in a convex form by using the Fourier series of the autocorrelation r % n :- 
Tjkgkgk-n of the filter g 13411 . Since r g; „ = r g _„ (real-valued symmetric sequence) we write 

L-l 

f g (v):=J]r g ,„0„(y) = |g(v)| 2 (17) 

n=0 

on the frequency band [— ^-] by using the basis (f> :— {I,2cos(27rv7o),2cos(27rv27o), . . . } and get |/3(v)| 2 = 
FgM ' I<?(k)| 2 - Due to the symmetry of </>„ and Sp we can restrict the constraints in < TT3T > to F — [0, and obtain the 



following semi-infinite linear problem: 

L-i 

max } r„„c„ such that < f„(v) < Al(v) for all ve 

s «=0 

Cn'= I \q(v)\ 2 </>n(v)dv. 



2Tn 



with M(v) := ^ F ^ y ^ MK \ ( . 



\q(v)f 



(18) 
(19) 



Since the FCC mask is piecewise constant, we separate M(v) into five sections vH(v) B35I1 and get the inequalities 



Vi = 1 5 : f g (v) < Mi(v) for v e [a it Pi\ 

with/3i = 1.61, p 2 = 1.99,/3 3 = 3.1,/3 4 = 10.6,/3 5 = 14 and ai = a 2 = ai = a A = 0,a 5 = /3 4 in GHz, seeFig.0 
The necessary lower bound for r g reads 



(20) 



f g (v) > for v e 



°' 2T 



= [0, 14]GHz. 



(21) 
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To formulate the constraints in ( fT9l ) for r g as a positive bounded cone in R L we approximate vH(v) by trigonometric 
polynomials ^ r,(v) := y,,„0„(v) of order L in the L 2 -norm iHl . The semi-infinite linear constraints in d20b describe 
a compact convex set [ " 

L-l 



34, (40),(41)]. To see this, let us introduce the following lower bound cones for 6 e [0, -^r] 



#iow(0) = \ r g € R L J] r g ,„0„(v) > 0, v € 



fe=0 



1 

27^ 



(22) 



For = the positive cone Kq = AT] ow (0) defines the lower bound in ( |2"TT i if we set To = 28GHz . To formulate the 
non-constant upper bounds, one can use the approximation functions F,(v) 13511 given in the same basis as |g(v)| 2 . 
For each i £ { 1 , . . . , 5} the bounds in (EOt are then equivalent to 

L 



2^(7i,n - r&n)<p n (y) > for v e [a,,/?,]. 



(23) 



For the upper bounds, we just have to set p, „ := y,- „ - r % n for each z = 1, . . . , 5 and n > 1, which leads to the upper 
bound cones 

L-l 



K uv (0i) = {r g e R L J]Pt»*»( v ) ^ °> v 6 



H=() 

L-l 



1 

9/, 

2T 



^upCft) = (r g s R L ^pi,„0„(v) > 0, v e [0, ft] 

n=0 

The five upper bound cones Kj are then 



Vi = 1 4 : := and ^ := ^ up (a 5 ). 



(24) 
(25) 

(26) 



Since the autocorrelation has to fulfill all these constraints, it has to be an element of the intersection. After this 
approximatiorj^ we g et from (fT3T > the problem 



max 



(27) 



n=0 



This is now a convex optimization problem of a linear functional over a convex set. By the positive real lemma [34J], 



these cone constraints can be equivalently described by semi-positive-definite matrix equalities, s.t. the problem d27l > 



3611 . The filter is obtained by a spectral factorization of 



is numerically solvable with the MATLAB toolbox SeDuMi |8, 
r g . Obviously g is not uniquely determined. 

Note that this optimization problem can also be seen as the maximization of a local L 2 -norm, given as the NESP value, 
under the constraints of local L°° -norms. 



4 Since the FCC mask divided by the Gaussian power spectrum is monotone increasing from to 10.6GHz we can let Ti, . . . , T4 overlap. 

5 The T, are approximations to the FCC mask with a certain error. Also, To is now fixed via the frequency range F. If one wants to reduce Tq, 

one has to reformulate the cones, hence y, and extend the frequency band constraints. Increasing To above 1 /28GHz is not possible, if one wants 
to respect the whole mask. 
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4. Orthogonalization of Pulse Translates 

In flQ 

a sequential pulse optimization was introduced, which produces mutually orthogonal pulses p° m = 
q *' Tg g„„ i.e. {p° iv p° n ) = S nm . Here each pulse p° n is generated by a different FIR filter g,„ e M L , which depends on the 
previously generated pulses p\, . . . , p° m _ x and produces pulses in L 2 ([-T q /2, (L-l)To+T q /2]). This approach is similar 
to the Gram-Schmidt construction in that it is order-dependent, since the first pulse p° can be optimally designed to 
the FCC mask without an orthogonalization constraint. We will now present a new order-independent method to 
generate from a fixed initial pulse p a set of orthogonal pulses \p° m }. Therefore we introduce a new time-shift T > 0, 
namely the PPM shift in (0, to generate a set of N - 2M + 1 translates {/?(■ - mT)}^ = M , i.e. M shifts in each time 
direction. The orthogonal pulses are then obtained by linear combinations of the translates of the initial pulse p. For a 
stable embedding of the finite construction we restrict the initial pulses to the set Lj := L 2 {[-T p /2, T p /2]) of centered 
pulses with finite duration T p . To study the convergence we need to introduce the concept of regular shift-invariant 
spaces. 

4.1. Shift-Invariant Spaces and Riesz-Bases 

To simplify notation we scale the time axis so that T = 1. Let us now consider the set So(p) := span(p(- — n)} of all 
finite linear combinations of \p{- - n)}, which is certainly a subset of L 2 . The L 2 -closure of So{p) is a shift-invariant 
closed subspace S(p) := So(p) c L 2 , i.e. for each / e S(p) also {/(■ - n)) c S(p). Since S(p) is generated by a single 
function p we call it a principal shift-invariant (PSI) space and p the generator for S{p). In fact, S{p) is the smallest 
PSI closed subspace of L 2 generated by p. Of course not every closed PSI space is of this form Jj^l . In this work we 
are interested in spaces which are closed under semi-discrete convolutions dl~2T > with ( 2 sequences, i.e. the space 

V(p):={p*'c \ cef 2 ) (28) 

9 



endowed with the L 2 - norm. Note that V(p) is in general not a subspace or even a closed subspace of L 2 Il37ll38il . But 
to guarantee stability of our filter design V(p) has to be closed, i.e. has to be a Hilbert subspace. More precisely, the 
translates of p have to form a Riesz basis. 



Definition 1. Let fi be a Hilbert space. \e„} c 'H is a Riesz basis for span{e„] if and only if there are constants 
< A < B < oo, s.t. 

2 



A llclli < 



<B||c||^ 2 for all c e I 2 . 



(29) 



■H 



In this case span{e„) becomes a Hilbert-subspace of U. For SI spaces in L 2 — "H we get the following result. 



Proposition 1 (Prop.l in 1 39]). Let p e L 



Then V(p) is a closed shift-invariant subspace ofL 2 if and only if 

A \\c\\ 2 e _ < \\p *' c\f L2 < B \\cf e for all cet 2 
holds for fixed constants < A < B < oo. Moreover, [p(- — n)) is a Riesz-basis for V(p). 



(30) 



If the generator p fulfills ( [30) . then V(p) - S{p) by [40] and we call p a stable generator and V(p) a regular 
PSI space [37]. An orthonormal generator {square-root Nyquist pulse, shift-orthonormal pulse) p° for V(p) is a 



generator with (p°(- - n),p°(- - m)) - 5 mn for all n, m e Z [ 1 1? ]. Benedetto and Li [41] showed that the stability and 
orthogonality of a generator p e L 2 can be described by the absolute [0,1] -integrable periodic function O p e L 1 ( [0, 1 ] ) 
of p defined for v almost everywhere (a.e.) as 



k 

They could show the following characterization HQ. 



(31) 



Theorem 1 (Th. 7.2.3 in H42I0 . A function p e L 2 is a stable generator for V(p) if and only if there exists < A < 
B < oo such that 



A <0,,(v) < B for v a.e. 
and is an orthonormal generator for V(p) if and only if 

®p( v ) — 1 for v a.e.. 



(32) 



(33) 



Proof. For a proof see Th. 7.2.3. (ii) and (iii) in 114211 . In our special case we have B = 1, Note, that the Riesz sequence 
and orthonormal sequence are bases for their closed span, meaning that in our case Sip) = V(p). □ 

Due to this characterization in frequency there is a simple "orthogonalization trick" for a stable generator given 
in (G3, which was found by Meyer, Mallat, Daubechies and others J^JLlQ. Prop. 7.3.9]. Unfortunately, this does 
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not provide an "a priori" construction in the time domain and does not lead to a support control of the orthonormal 
generator in time, as necessary for UWB-IR. 

Contrary to an approximation in the frequency-domain we approach an approximation in time-domain via the 
Lowdin transformation. We will show that in the limit the Lowdin transformation for shift-sequences is in fact 
given in frequency by the orthogonalization trick (Q~|i. By using finite section methods we establish an approximation 
method in terms of the discrete Fourier transform (DFT) to allow an easy computation. Furthermore, we show that the 
Lowdin construction for stable generators is unique and optimal in the L 2 -distance among all orthonormal generators 
and corresponds to the canonical construction of so called tight frames (given later). 

4.2. Lowdin Orthogonalization for Finite Dimensions 

Since the Lowdin transformation is originally defined for a finite set of linearly independent elements in a Hilbert 
space < H, we will use the finite section method to derive a stable approximation to the infinite case. For this we consider 
for any M e N the symmetric orthogonal projection Pm from £ 2 to l 2 M — {c e l 2 \ supp c c {-M, . . ., M}} defined for 
c e t 2 by PmC := c M = (0, . . . , 0, c_m, • ■ ■ , cm, 0, . . . , 0). Then the finite section Gm of the infinite dimensional Gram 
matrix G of p e L 2 , given by 

[G]„ m := <p(- - m), p(- - n)> for n, m e Z, (34) 

can be defined as Gm := PwGPm, see J44I Prop. 5.1.5]. If p satisfies ( f29T > and if we restrict the semi-discrete 
convolution p *' c to fi M , we obtain a 2M + 1 dimensional Hilbert subspace V M {p) of V{p) = < H. Then the unique 
linear operation B M , which generates from {p(- - m)}^ = M an orthonormal basis (ONB) {Pm M }„ = _ M for V M (p) and 
simultaneously minimizes 

M 

2 \\B M p(--m)-p(--m)\\ 2 L2 (35) 

m=-M 

is given by the (symmetric) Lowdin transformation 



and can be represented in matrix form as 



M 

plf := B M p(- -m)=Yj [G* WO " ») for a11 me{-M,..., M], (36) 

n=-M 

where we call each p°^ M a Lowdin orthogonal (LO) pulse or Lowdin pulse. 

1 _ 1 

Here G M 2 denotes the (canonical) inverse square-root (restricted to ( 2 M ) of Gm- Note that G M 2 is not equal to 

PmG ^Pm- Since the sum in <f36b is finite, the definition of the Lowdin pulses is also pointwise well-defined. In the 

next section, we will see that this is a priori not true for the infinite case. If we identify the corresponding m'th row 

of the inverse square-root of Gm with vectors bjjf = ([G m 2 ],„_m, • ■ • , [G m 2 ]„ + m) we can describe ( 1361 ) by a FIR filter 

bank as 

PT = P*'K . me\-M,...,M). (37) 
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Unfortunately, none of these Lowdin pulses is a shift-orthonormal pulse, which would be necessary for an OOPPM 
transmission. In the next section we will thus show that the Lowdin orthogonalization converges for M to infinity to 
an IIR filter b given as the centered row of G " . This IIR filter generates then a shift-orthonormal pulse, namely the 
square-root Nyquist pulse defined in ([]]). Hence the Lowdin orthogonalization (l36l l provides an approximation to our 
OOPPM design. In the following we will investigate its stability, i.e. its convergence property. 

4.3. Stability and Approximation 

In this section we investigate the limit of the Lowdin orthogonalization in (l36l l for translates (time-shifts) of the 
optimized pulse p with time duration T p < oo where we further assume that p is a bounded stable generator. If we set 
K := [T p ] then certainly p e L 2 K . In this case the auto-correlation of p 



r p (t) := (p * p_)(t) = J^ p (t)p(t - f)dr, 



(38) 



with the time reversal p-(t) := p(— t) is a compactly supported bounded function on [-K, K]. Due to the Poisson 
summation formula we can represent <& p almost everywhere by the Fourier series ( TToT l (7q = 1) of the samples \r p {n)} 

oo K 

<W = J] r p (n)e- 2 " inv = £ [GUe^, (39) 

n=— oo n=-K 

which is the symbol of the Toeplitz matrix G, since we have from (l34l i and d38l that r p (n - m) = [G]„,„. Moreover the 
symbol is continuous since the sum is finite due to the compactness of r p . 

On the other hand the initial pulse p is a Wiener function^ J44I Def. 6.1.1] so that O p defines a continuous 
function and condition d32l > holds pointwise J47I Prop.l]. Since both sides in d39i > are identical a.e. they are identical 



everywhere by continuity (see also [44, p. 105]). Thus, the spectrum of G is continuous, strictly positive and bounded 



by the Riesz bounds. Hence the inverse square-roots of G and G« exists s.t. for any MeN (by Cauchy's interlace 



theorem, [48, Th. 9.19]) 



A\ M <G M <B\ M and il v < G M ' < \l„, (40) 

B A 

where 1m denotes the identity on ( 2 M . Now we can approximate the Gram matrix by Strang's circulant precon- 



ditioner [49], s.t. the diagonalization is given by a discrete Fourier transform (DFT) 115011 . To get a continuous 
formulation of the approximated Lowdin pulses we use the Zak transform 05111 . given for a continuous function / as 

(Zmy)-=Yjf(t-n)e Mnv for f,veR. (41) 
Our main result is the following theorem. 



1 Wiener functions are locally bounded in L°° and globally in I . 
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Theorem 2. Let K € N and p e L 2 K be a continuous stable generator for V(p). Then we can approximate the limit set 
of the Lowdin pulses {p° m } by a sequence of finite function sets {p° 1 ] M } ! ^ l= _ M , which are approximate Lowdin orthogonal 
(ALO). The functions p„ M are given pointwise for M > K and m e [-M, . . . , M) by the Zak transform as 

f —Iniml , 



such that for each m € Z 







P°Jf) = lim ~ P :f{t) 

M— >oo 



l<M-# 



(42) 



else 



(43) 



converges pointwise on compact sets. The limit in ( I43l l can foe stated as 

p°(v) = p(y) • (0 p (y)) 4 (44) 
for v e M in the frequency-domain. Hence the Lowdin generator p° :— p° G is an orthonormal generator for V(p). 

Proof. The proof consists of two parts. In the first part we derive an straightforward finite construction in the time 
domain to obtain time-limited pulses d42l i being approximations to the Lowdin pulses. Using Strang's circulant 
preconditioner the ALO pulses can be easily derived in terms of DFTs. In the second part we will then show that 
this finite construction is indeed a stable approximation to the square-root Nyquist pulse. Here we need pointwise 
convergence, i.e. convergence in (the set of bounded sequences). Finally, to establish the shift-orthogonality we 
use properties of the Zak transform. 

Since the inverse square-root of aN x N Toeplitz matrix is hard to compute, we approximate for any M > K the 



Gram matrix Gm by using Strang's circulant preconditioner Gm [49 



and banded such that we can define the elements of the first row by [53 



5211 . Moreover, the Gram matrix is hermitian 



(4.19)] as 



[G 



AfJOn •- 



r p (n) 



n€{0,...,K] 



rJN-n) ne{N- K,...,N- 1} 



(45) 







else 



(46) 



Here we abbreviate := 2M + 1. The crucial property of Strang's preconditioner Gm is the fact that the eigenvalues 
A/(Gm) are sample values of the symbol fl> p in OTb . This special property is in general not valid for other circulant 
preconditioners |52[. To see this, we derive the eigenvalues by J53I Theorem 7] as 

K N-l 

;if :=^GM) = 2 r '' (n)e ~ 27 "' 4 + X r pW~n)e- Mn * for I 6 {0, . . . ,2M] 

«=0 n=N-K 

by inserting the first row of Gm given in d45l ). If we set in the second sum n' — n - N we get from d39b 

n=-K 
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If 



2M+1 



(47) 



Since p is compactly supported the symbol <t> p is continuous and the second equality in ( l47l i holds pointwise. More- 
over, the Riesz bounds d32l of p guarantee that Gm is strictly positive and invertible for any M. Now we are able to 
define the ALO pulses in matrix notation by setting in p M (0 := {/?( f_ »))f _ M = (p( t+M )> • • • > P(t~M)) T for any t e K 



p°' M (f) := G M I p M (f) = F M D M 5 F^p M (f), (48) 
since the circulant matrix G« = F^D^F*^ can be written by the unitary N x N DFT matrix Fm, with 



[F M ]„„, := -^— e - 2ni Tr with n,m€{0, . . . , 2M) 



(49) 



and the diagonal matrix Dm of the eigenvalues of Gm- 

Let us start in ( l48b from the right by applying the IDFT matrix F* M , then we get for any kth component with 
k e {0, . . . , 2M) 

M 2M 

(50) 



[F M p M (f)]i 



1 / \ 

—\Y j e 7 * i i k -p(t + M-n)+ Yj e 2ni ^ k ■ p{t + M - n)\ 



r i 



A/ 



2/n ■ ■ 



Pit -J)- 



(51) 



Next we multiply with the components [D M 2 = 5^// -yj/l, of the inverse square-root of the diagonal matrix Dm 

{ 



D,/F M p M (r) 



=— r 



„2m L -r r k 



Pit -J) 



1 



e^'-pit-j) 



i=-M 



(52) 



In the last step we evaluate the DFT at m e {0, . . . , 2M) 

-2mm- A Jl 



2 A/ 

z 

/=() 



2> 



iV 



2A/ 

z> 

/=() 



-2ml- 



XjL-M Pi { ~ j)e " 



(53) 



where we used the Zak transform (PTTT i of r p to express the eigenvalues if. In the next step we extend the DFT sum 
of the numerator in d53l to an infinite sum. This is possible since p(- - k) always has the same support length K for 
each k e Z. Thus, for all \t\ > M - § the non-zero sample values are shifted in the kern of Pm; hence p M (f) = 0. On 
the other hand for \t\ < M - j any shift \j\ > M results in p(t - j) = 0. If we also set for each M the index k := m — M 
in 031 then the continuous ALO pulses can be written as 



Pk W 







< M-# 



(54) 



else 



7 By slight abuse of our notation we understand in this section any matrix as an iV X iV matrix and p M (0, P°' M W as TV-dim. vectors for any (El. 
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This defines for each k the operation Bf by p° k ' M = B k p. For each t e R the term y^KO^) is a continuous function 
in v since the Zak transforms are finite sums of continuous functions and the nominator vanishes nowhere. This is 
guaranteed by the positivity and continuity of <J) p due to (|40T >. Hence the ALO pulses are continuous as well. 

The second part of the proof shows the convergence of our finite construction to a square-root Nyquist pulse. 
Therefore we use the finite section method for the Gram matrix. Gray showed J53I Lemma 7] that 



(Zp)(t,v) 



Gm - Gj 







(55) 



as M — > oo in the weak norm ||A||^, :- 1/N Yjj tf(A) implying weak convergence of the operators. Since Gm is strictly 
positive for each MeNwe get by [54] 



0. 



Unfortunately this does not provides a strong convergence, which is necessary to state convergence in I 2 

G m PmC— >G 5 c foranyceT. 



(56) 



(57) 



However, from [55] finite strong convergence can be ensured, i.e. convergence of ( 157b for all c € t 2 M , for each 
M' e N. But for any ( e 1 there exists an M' sufficiently large, due to the compact support property of p, such 
that c = p(f) := {p(t - n)} e t 2 M ,. This is in fact sufficient, since it implies pointwise convergence in d^., of (IBTt . i.e. 
component- wise convergence for each t e R. Let us take for each t e Rthe number M' e N such that max] \t\ ,K} < M' . 
Then we can define the limit of the k\h component as 



pl(t) := lim p° k M {t). 

M— >oo 



If we define Av 



l 



, i; j and v; = /Av we can write for d58l l by inserting ( f54t 

2M „-2mkv, ( 



lim pl' M (t) = lim V 



e- Zmkv >(Zp)(t,vi) 



Av. 



j_ ^(Zr p )(0,v,) 

n 

Using the quasi-periodicity 115 1L (2. 18), (2. 19)] of the Zak transform for leZwe have for any f, v e 

(Zp(- - k))(t, v) = (Zp)(t -k,v) = e- 2 " ivk {Zp){t, v). 



We can express the partial sum on the right hand side of d59l ) in the limit as a Riemann integral for each t e 



p° k (t) := lim p° k '"(t) = f 
«^ Jo 



„ aMfA _ , (Zp{--k))(t,v) 



dv. 



(58) 



(59) 



(60) 



(61) 



V(Zr p )(0,v) 

This shows that Bfp converge pointwise for M — > oo to B^p = Bop(- - k) — p° k for each k. The sequence {p° k } is 
then generated by shifts of the centered pulse p° := p° since the shift operation commutes with B := Bq. This in turn 
commutes with the Zak transformation, i.e. for all t e R we have 



(Bp)(t-k) 



I 



(Zp)(t-k,v) 
V(Zr„)(0,v)" 

= (Bp(-k))(t). 
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dv : 



(Zp(-kW,v) 
V(Zr p )(0,v) 



dv 



(62) 
(63) 



From ( I6TT ) it is now easy to show that p° is an orthonormal generator. We write the left hand side of ( |6*TI ) in the Zak 

3lfl t 



domain, by applying the Zak transformatior|j to p° 



(Zp )(f,v) = :. (64) 

V(Zr p )(0,v) 



(65) 



If we multiple ( l64l l by the exponential and integrate over the time we yield for every v e K 

Jo Jo JZrp(0, v) 

Since <D p = (Zr p ){0, •) is time-independent we get the "orthogonalization trick" (fTJ by using in d65l l the inversion 
formula 115 lL (2.30)] of the Zak transform 



P°(y) =P(y) • (a>„(v))~ a = p°(v). (66) 

Again, this is also defined pointwise since the right hand side is continuous in v. It can now be easily verified that p° 
fulfills the shift-orthonormal condition (1331 . which shows that p° is an orthonormal generator for V(p). □ 



Remark. Note, that relation ( f60b induces a time-shift. To apply this to the ALO pulses in (l54l the time domain has to 
be restricted further. Hence the ALO pulses do not have global shift character for finite M e N, but locally, i.e. p° k ' M 
shifted back to the center matches p° M for t e [-M + § + |Jt| , M - f - |jfc|]: 

~°* ( t + lA 1 V -2mU ^n-MP(t + k-n)e^ n I v Z™_ M p(t + k - ri)e^ L ^ 

ft» + «- s 2> ^ - S Z ^ <«' 

— # — ■ 

Since p(r + M + |Jk|) = and p(t - M - \k\) = for \k\ < M and \t\ < M - f - |Jk|, we end up with 

^* ( f + ^) = 77/, p= =P W- (69) 

^ ' ft 

For all < M the ALO pulses have the same shape in the window \t\ < M — \ — if we shift them back to the origin. 

Moreover, the ALO pulses are all continuous on the real line, since they are a finite sum of continuous functions 
by definition d48t . Hence each ALO pulse goes continuously to zero at the support boundaries. So far it is not clear 
whenever p° is continuous or not. Nevertheless its spectrum p° is continuous and so we can state p° = p° almost 
everywhere. Hence the orthogonalization trick defines the square-root Nyquist pulse p° only in an L 2 sense. 



A similar result is also known in the context of Gabor frames, see also M, 8.3] 



so Q 
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5. Discussion of the Analysis 

In this section we will discuss now the properties of our OOPPM design for UWB, i.e. the optimization and 
orthogonalization, which can be completely described by an IIR filtering process. First we will relate the Lowdin 
orthogonalization to the canonical tight frame construction. Afterwards we will show in Section [5721 that the Lowdin 
transform yields the orthogonal generator with the minimal L 2 -difference to the initial optimized pulse. This is the 
same optimality property as for canonical tight frames |20]. But such an energy optimality does not guarantee FCC 
compliance. So we will discuss in Section IB31 the influence of a perfect orthogonalization to the FCC optimization. 
Finally, we will discuss the implementation of a perfect orthogonalization by FIR filtering. 

5.7. Relation Between Tight Frames and ONBs 

Any Riesz basis [pk\ for a Hilbert space fi is also a exact frame for 'H with the frame operator S defined by 

S:<H^<H, f»Sf= (70) 

k 

where the frame bounds are given by the Riesz bounds 0<A < B<oo of [pk] [42, Th. 5.4.1, 6.1.1], i.e. 

A ll/II* <<£/,/>„ <fi||/ll* for any / e "K. (71) 
Here (•, -), H denotes the inner product in < H and ||-||^ the induced norm. Since S is bounded and invertible, i.e. the 



inverse operator exists and is bounded B42I1 . we can write each /eWas 

f = SS- 1 f = Yj,S- 1 f,Pk) n Pk. (72) 

k 

In this case the Lowdin orthonormalization corresponds to the canonical tight frame construction. 



Lemma 1. Let the sequence {p^} be a Riesz basis for the Hilbert space 'H := spanj/?^} and G its Gram matrix. Then 
the canonical tight frame {p° k ) is given for each k e Z by: 

p° k :=S-hk = J][G-huPi (73) 

in an L 2 -senseQ 



Proof. See Appendix A 



If we now set p„ :- p(- - n) e L 2 the Riesz-basis is generated by shifts of a stable generator and "H = V(p) 
becomes a principal shift-invariant (PSI) space, 

which is a separable Hilbert subspace of L? as discussed in Section |4~T~I The canonical tight frame construction 
then generates a shift-orthonormal basis, i.e. an orthonormal generator. The reason is that shift-invariant frames and 



9 This statement was already given without further explanation by Y. Meyer in [43] equation (3.3). Note that Y. Meyer used condition (3.1) and 
(3.2) in 1 43] which are equivalent to the Riesz basis condition. 
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Riesz bases are the same in regular shift-invariant spaces 1 56, Th.2.4]. So any frame becomes a Riesz basis (exact 
frame) and any tight frame an ONB (exact tight frame). Hence for regular PSI spaces there exists no redundancy for 
frames. This generalize the Lowdin transform for generating a square-root Nyquist pulse to any stable generator p. 

From Meyer 14311 we know that ( 1731 ) can be written in frequency domain as the orthogonalization trick. Therefore 
the limit of the Lowdin transformation B - B: V(p) — > V(p) 



f^Bf 



f 

Jo 



(Z/)(-,v) 
VZi»(0,v)" 



■xlv 



(74) 



equals the inverse square-root of the frame operator in ( f73l . 



5.2. Optimality of the Lowdin Orthogonalization 

Janssen and Strohmer have shown in 112011 that the canonical tight-frame construction of Gabor frames for L 2 is 
via Ron-Shen duality equivalent to an ONB construction on the adjoint time-frequency lattice. Furthermore they have 
shown that among all tight Gabor frames, the canonical construction yields this particular generator with minimal 
L 2 -distance to the original one. However, for SI spaces this optimality of the Lowdin orthogonalization has to be 
proved otherwise. To prove this we use the structure of regular PSI spaces. 

Theorem 3. The unique orthonormal generator with the minimal L 2 distance to the normalized stable generator 
p G L 2 for V(p) is given by the Lowdin generator p° . 



Proof. See | Appendix B| 

Nevertheless, We have to rescale the orthonormal generator p° to respect the FCC mask, see Section|6] For this 
the maximal difference of the power spectrurrQ of the (normalized) optimal designed pulse and the orthonormalized 
pulse is of interest, i.e. 



'| 2 - l^°| 2 || L « = ess sup 



<D»(v) - 1 



4>p(v) 



|p(v)|' 



3>„ 



4>„ 



(75) 



This shows again that this L°° distortion is also determined by the spectral properties of the optimal designed pulse p 
and its Riesz bounds. Unfortunately it is very hard to control the optimization and orthogonalization filter simultane- 
ously as will be shown in the next section. 



5.3. Interdependence of Orthogonalization and Optimization 

The causal FIR operation in ([T3l of a fixed initial pulse q of odd order L with clock rate 1/7*0 can be also written 
in the time-symmetric form as a real semi-discrete convolution 



P = ?*'r„g for ge^L(K). 



(76) 



In fact the L™-distance of the FCC mask Sp and in F is relevant, assumed \p°\ is bounded by Sp. 
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In this section we investigate the interdependence of the IIR filter b and the FIR filter g, i.e. the interdependence of 
the orthogonalization filter in Section l4~2l and the FCC optimization filter in Section|3]for different clock rates. So far 
we have first optimized spectrally and afterwards performed the orthonormalization. In this order for a chosen q, the 
orthogonalization filter b depends on g, hence we write b = b ? . Moreover the clock rates of the filters differ, hence 
we stick the time-shifts as index in the semi-discrete convolutions. For the r-shift-orthogonal pulse we get then 



p T '° = (q *' To g) 4 b g 



(77) 



Let us set T = A7o and T q = N q To for N q e N, A > 0. Since the filter clock rate of g is fixed to l/7o to ensure 
full FCC-range control, the variation is expressed in A. To get rid of To we scale the time t to ?' = f/7o such that 
the time-shift of g is 7^ = 1. We observe the following effects: 1) If i e N then p A -° = q A -°. 2) If A e N then the 
distortion by o g is limited periodically to the interval [-jt, sj]. 3) If q is already 1 -shift-orthogonal and A e N then 
b ? can be omitted and instead adding an extra condition on g to be a A-shift orthogonal-filter, i.e. r g {kA) = dto, which 
ensures A-shift orthogonality of the output p. To see point 1), let us first orthogonalize q by b 9 and ask for the filter 
g A „ which preserves the A-orthogonalization in the presence of g. Hence we aim at 



»A,o 



P ' =P *A b g = 1 

But from (l44l > we know how b^ acts in the frequency-domain: 



A.o 



»A,o 



(v)| 2 = 



\P(vW 



- I \q & '°(v) ■ g A , g (v)| ■ 



(78) 



(79) 



Since |p(v)| 2 = |g(v) • q(v)\ 2 and jeN we get by the ^-periodicity of |g|' 



|§ A '°(V) • § A , g (v)| 2 



lg(v)| 2 • \q(v)\ 2 



\q(y)\ 2 



7 A '°(v)f 



(80) 



ilg(v)l 2 -^|§(v-|)| {lk\q(y-j)\ 
Hence we get § A ,g(^) = £*o> which shows 1). The price of the orthogonalization is the loss of a frequency control, 
since the frequency property is now completely given by the basic pulse q and time-shift A. In Fig. [10] the effect is 
plotted for A e [1,2] and L — 25. For small A the distortion is increase by the orthogonalization. This also shows that 
a perfect orthogonalization and optimization with the same clock rates is not possible. 

In 2) a perfect orthogonalization does not completely undo the optimization, since T — A > Tq — 1 . For A = 2 we 



can describe the filter by using the addition theorem in |g(v + 1 /2)| = f g (v + 1/2) = 2r gj o - ?g(v) by 



z 



P|v-- 



f„(v) 



z 



q v + 



2k 



f„(v) 



Z 



q v + 



2Ar+l 



? (v)[0> 9 (v) 



2r, 



f„(v) 



1 W(y)] (81) 



=:0'(v) 



which results in the filter power spectrum ( |79l ) 

|g2, g (v)| 2 =(l+( 



2r, 



g,o 



f g (v) 
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- 1 



a> a (v) 



(82) 



But since we fixed A = 2 and q we can calculate <b q ,<&' and q A, °. This provides a separation of the filter power 
spectrum f g (v) = |g(v)| 2 and the orthogonalization. Unfortunately, this does not yield linear constraints for r g . 
Finally, note that the time-shifts and hence the filter clock rates have to be chosen such that an overlap of the basic 
pulse occurs. Otherwise a frequency shaping is not possible. 

Case 3) assumes already a shift-orthogonality. We only have to ensure that the spectral optimization filter g 
preserves the orthogonality. This results in an extra orthogonality constrain for the filter g, which can be easily 
incorporated in the SDP problem of Section [3] see Davidson et. al in [? ]. 

Summarizing, the discussion above shows that joint optimization and orthogonalization is a complicated problem 
and only in specific situations a closed-form solution seems to be possible. 

5.4. Compactly Supported Orthogonal Generators 

For PPM transmission a time-limited shift-orthogonal pulse is necessary to guarantee an ISI free modulation in a 
finite time. Such a PPM pulse is a compactly supported orthogonal (CSO) generator. In PPM this is simply realized 
by avoiding the overlap of translates. 

To apply our OOPPM design it is hence necessary to guarantee a compact support of the Lowdin generator p° 
given in Theorem|2] In this section we will therefore investigate the support properties of orthogonal generators. The 
existence of a CSO generator (with overlap) was already shown by Daubechies in Jl^l . Unfortunately, she could not 
derive a closed-form for such an CSO generator. Moreover, to obtain a realizable construction of a CSO generator 
this construction has to be performed in a finite time. So our Lowdin construction should be obtained by a FIR filter. 

PSI spaces of compactly supported (CS) generators, were characterized in detail by de Boor et al. in and 
called local PSI spaces. If the generator is also stable, as in Theorem [2] then there exists a sequence c e I 1 such 
that p *' c is an CSO generator. Moreover, any CSO generator is of this form. To investigate compactness, de Boor 
et.al. introduced the concept of linear independent shifts for CS generators. The linear independence property of a 
CS generator p is equivalent by J37I Res. 2.24] to 

{(Jf/>)(z-n)}„EZ*0 for all zeC (83) 

where (Jzfp) denotes the Laplace transform of p. This means (J£p) do not have periodic zero points. Note that this 
definition of independence is stronger than finitely independence, see definition in If we additionally demand 
linear independence of p in our Theorem [2] then this CS generator is unique up to shifts and scalar multiplies. 
Furthermore, a negative result is shown in [37], which excludes the existence of a CSO generator if p itself is not 
already orthogonal. But if p is already orthogonal, then p is unique up to shifts and scalar multiplies and then the 
Lowdin construction becomes a scaled identity (normalizing of p). The statement is the following: 

Theorem 4 (Th. 2.29 in [37]). Let p e L 2 be a linear independent generator for S{p) which is not orthogonal, then 
there does not exists a compactly supported orthogonal generator p° for S(p), i.e. there exists no filter b € t 1 such 
that p° = p *' b. 
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If p is a linear independent generator which is not orthogonal, then the Lowdin generator p° has not compact 
support. We extend this together with the existence and uniqueness of a linear independent generator for a local PSI 
space S(p): 

Corrolary 1. Let p e L 2 be compactly supported. If there exists a compactly supported orthogonal generator p° for 
S(p), then it is unique. 



Proof. Any CSO generator p° e L 2 is a linear independent generator by 1371 Prop. 2.25(c)]. Since the linear indepen- 



Remark. In any case there exists an orthogonal generator for a local PSI space. For a stable CS generator p our 
Theorem [2] gives an explicit construction and approximation for an orthogonal generator by an IIR filtering of p. If 
the Lowdin generator is not CS, it is the unique orthogonal generator with the minimal L 2 -distance to the original 
stable CS generator by Theorem[3] So far it is not clear whether there exists a IIR filter eel 2 which generates a CSO 
generator from a stable CS generator or not. What we can say is that if the inverse square-root of the Gram matrix is 
banded, then the rows corresponds to FIR filters which produce CSO generators, since the semi-discrete convolution 
reduces to a finite linear combination of CS generators. So this is a sufficient condition for the Lowdin generator to 
be CSO, but not a necessary one. 

6. Application in UWB Impulse Radio Systems 

Here we give some exemplary applications of our filter designs developed in Section[3]and Section|4]for UWB-IR. 

FIR filter realized by a distributed transversal filter. The FIR filter is completely realized in an analog fashion. It 
consists of time-delay lines and multiplication of the input with the filter constants. Note also that these filter values 
are real- valued. An application in UWB radios was already considered in I57L 

Transmitter and receiver designs. Our channel model is an AWGN channel, i.e. the received signal r(t) is the trans- 
mitted UWB signal u(f) given in © by adding white Gaussian noise. For simplicity of the discussion we omitted the 
time-hopping sequence (c„) in (0. We propose now three N - ary waveform modulations for our pulse design. Since 
our proposed modulations are linear and performed in the baseband, the signals (pulses) are real-valuedl^ 

(a) A pulse shape modulation (PSM) with the Lowdin pulses iri'°' M )f = _ J ,, which corresponds to a A^-ary orthogonal 
waveform modulation. The n'th message mis transmitted as the signal u(t) = yf& a n p T ^' M (t—nT s ) - a n s m {t—nT s ). 
The receiver is realized by N correlators using the Lowdin pulses as templates. From the correlators output l m the 
absolute value is taken due to the random amplitude flip by a„, see Fig. [3] 




dent generator is unique by 1371 Th. 2.28(b)], the CSO generator is as well. 



□ 



Note that our proposed pulse design can be also used for a complex modulations (carrier based modulation), e.g. for OFDM or FSK. 
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Figure 3: Optimal Receiver for JV-ary orthogonal modulation with scaled Lowdin pulses .v,„ := ~j£^ l pJ l f' M . 
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Figure 4: Matched filter receiver for an ALO and LO-OPPM scheme. 

(b) The centered ALO and LO pulse p T >°> M := pl'° M resp. p T -°> M := M for an OPPM design are in fact a 
non-orthogonal modulation scheme with a matched filter at the receiver, see Fig. |4] The n'th message m is 
then transmitted as u{f) = V& a„p T -°' M (t - nT s - mT) in the ALO-OPPM design. The matched filter output 
37(f) = £l' +Ts/2 r(T)p T '°' M (T - t)dr is sampled for the n'th message at y(nT s + mT). (for LO-OPPM use p T -°- M ) 

(c) The limiting OOPPM design with the Lowdin pulse p° is not practically feasible, since we have to use an IIR filter. 
Hence we only refer to this setup as the theoretical limit. The transmitted signal would be u(t) = V£ a„p T '°(t - 
nT s - mT) with the matched filter h(f) — p T '°(t). Note that the receiver in Fig. [4] also has to integrate over the 
whole time due to the unlimited support, which would produce an infinity delay in the decoding process. 



Scaling with respect to the FCC mask. The operations B M and B M generate pulses which are normalized in energy 
but do not respect anymore the FCC mask. So we have to find for the mth pulse its maximal scaling factor a„, > s.t. 

|2 



\<X m -p°„ M (v)\~ <S F (V) 

is still valid for any v € F. This problem is solved by 



(84) 
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For the scaled Lowdin pulses we can easily obtain the following upper bound for the NESP value ( fT4t 



I, 



\iy* n°' M \ ll-°.M|l 2 
\ a mPm | \\p m , 

< - — 
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L~([0,F] 
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L~([0,F] 



1 



-vo,M|- 



(86) 



L°°qo,F] 



where we denoted with &p p — J F Sf the allowed energy of the FCC mask in the passband F p . Thus, the maximization 
of the symbol energy & under the FCC mask is the maximization of or* in d85l l. i.e. a maximization of the L°°-norm 
in the frequency domain. 

Remark. We want to emphasize at this section, that the FCC spectral optimization is rather an optimization of the 
power in an allowable mask given by the FCC, than an optimization of the spectral efficiency of the signal. Although 
the PSWF's have the best energy concentration in [-W, W] among all time-limited finite energy signals fistl . they are 
not achieving the best possible nesp value 02D. The reason is, that spectral efficiency with respect to UWB is to utilize 
as much power from the UWB window, framed by the FCC mask Sf, as possible. So the price of power utilization in 
[-W, W] under the UWB peak power limit Sf, is a loss of energy concentration and hence a loss of spectral efficiency 
compared to the PSWF's. 



6.1. Performance of the Proposed Designs 

For a given transmission design, consisting of a modulation scheme and a receiver, the average bit error probability 
P e over Eh /No is usually considered as the performance criterion. We consider real-valued signals in the baseband with 
finite symbol duration T s . The optimal receiver for a non-orthogonal A^-ary waveform transmission is the correlation 
receiver with M correlators, see Fig. |3]with maximum likelihood decision. 

N-ary orthogonal PSM for scheme © above. The average (symbol) error probability for N-ary orthogonal pulses 
with equal energy & can be upper bounded by 15811 



P e <(N- l)erfc 




(87) 



Note, that this error probability is the same as for an orthogonal PPM modulation (fTTb . To obtain equal energy symbols 
and FCC compliance we have to scale each Lowdin pulse with Vfi - a* - min m {ar^,}. 

N-ary overlapping PPM for scheme (JbJ above. Here we can substitute the N correlations by one matched filter 
h = p T '°' M resp. h - p T '°' M and obtain the statistics \l m \ by sampling the output. The average error probability per 
symbol P e given exactly in |5^, Prob. 4.2.1 1] for equal energy signals and can be computed numerically. The energy 
is given by Vfi = a* — a* and VI) = 5rS calculated in (f85T > for p T <°< M resp. p T '°' M . Upper bounds obtained in I59I can 
be used for the ALO resp. LO average error probability 



Pe< jjerfc 



1 N 

and P e < — ^ erfc 

2 j=2 



(88) 
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Figure 5: Optimal pulse p for L = 25 and basic pulse q = p for L= 1 in time-domain in generic units (GU). 



with p\j = &r p T.o,M(jT) and p\j = &r p T,°,M(jT), since the symbols are given by Sj = yf&p T '°' M (- - jT) resp. Sj = 
y[Ep T '°' M (- - jT) for j = l,...,N. The error probabilities depend on the pulse energy and on the decay of the sampled 
auto-correlation defined in ( f38l >. 

6.2. Simulation Results 

The most common basic pulse for an UWB-IR transmission is the Gaussian monocycle: q(t) t ■ exp(-f 2 /cr 2 ) 
where <x is chosen such that the maximum of \q( f)\ 2 is reached at the center frequency f c = 6.85GHz of the passband 
|4J]. Since we need compact support and continuity for our construction, we mask q with a unit triangle window A 
instead of a simple truncation. Also any other continuous window function which goes continuously to zero (e.g. the 
Hann window) can be used, as long as the lower Riesz bound A > can be ensured, see Theorem[2] We have used an 



algorithm in |3| to compute A and fi numerically. Note that for any continuous compactly supported function we ha\ c 
a finite upper Riesz bound B, see B60l Th.2.1]. The width (window length) is chosen to T q — — N q To as 0.21428ns, 



such that at least 99.99% of the energy of q is contained in the window [—T q /2, T q /2], see Fig. [5] We express all 
time instants as integer multiples of 7q. Also, in Fig. [5] we plot the optimal pulse obtained by a FIR filter of order 
L = 25 which results in a time-duration T p = 307o = 5T q of p. In our simulation we choose T q :- 6 To = N q To 
and L — 25 as the filter order of the FCC-optimization. Hence, the optimized pulses have a total time duration of 
T p = (L - l)7o + T q = 307o = N P T . 

The Riesz condition (l32l l has been already verified in [3] for this particular setup. Theorem[2]uses the normaliza- 
tion T = 1, Translating between different support lengths T' p — K is done by setting t := t'T p /K. Now the support 
of p(t') is [-K/2,K/2] with fixed T — 1. To obtain good shift-orthogonality, we have to choose M > K. This we 
control with an integer multiple m = 2, i.e. M = mK = 2K. The support length T s of all the LO (ALO) pulses is then 
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given as 



= (N - 1)7 + T„ = (2mK)T p IK + N p T 



(89) 



= 7>,.« = (2m + l)N p T = 150r . 



Now the time slot [-T p t,„.m /2, T p T,o,u/2] exactly contains N mutually orthogonal pulses {pm°' M }, i.e. N orthogonal 



symbols with symbol duration T s = T p t.°m having all the same energy and respecting strictly the FCC mask. This is 
a N-ary orthogonal signal design, which requires high complexity at receiver and transmitter, since we need a filter 
bank of N different filters. 

Our proposal goes one step further. If we only consider one filter, which generates at the output the centered Lowdin 
orthonormal pulse p T '°- M , we can use this as a approximated square-root Nyquist pulse with a PPM shift of T = T p /K 
to enable N-ary OPPM transmission by obtaining almost orthogonality. 

Advantages of the proposed design are: a low complexity at transmitter and receiver, a combining of g and b 
into a single filter operating with clock rate I/Tq and q as input if T — T p /{TqK) e N , a signal processing "On 
the fly" and finally a much higher bit rate compared to a binary-PPM. The only precondition for all this, is a perfect 
synchronisation between transmitter and receiver. In fact, we have to sample equidistantly at rate of 1 /T. The output 
of the matched filter h{—t) = p T '°' M (t) is given by 



and recovers the rath symbol. The statistic l m = y(mT) is the correlation of the received signal with the symbol s m , 
see Fig. |4] Note that the shifts have support in the window [—1.5T p t,„m, 1.5T p t.o,m], but are almost orthogonal outside 
the symbol window [-T p t.°,m /2, T p t,v,m/2] due to the compactness and approximate shift-orthogonal character of the 
Lowdin pulse p T, ° ,M . 

In Fig. [6] the centered orthogonal pulses for T = 5To = 5T q /6 match almost everywhere the original masked 
Gaussian monocycle, since the translates are almost non-overlapping, hence they are already almost orthogonal. For 
T = To — T q /6 the overlap results in a distortion of the centered orthogonal pulses, where the ALO pulses have high 
energy concentration at the boundary (circulant extension of the Gram Matrix). 

In Fig. [7] the pulse shapes in time for the centered Lowdin orthogonal pulses are plotted. The ALO pulses are 

matching the LO pulses for T > 2Tq almost perfectly such that we did not plot them, since the resolution of the plot is 

to small to see any mismatch. Only for the critical shift T — To a visible distortion is obtained at the boundary. In the 

next Fig. [8] we plotted therefore the ALO and LO pulse for T = 1.5Tq to show that the ALO pulses indeed converge 

very fast to the LO pulse if T » To- The reason is that for small time shifts of p the Riesz bounds and so the clustering 

_ i _i i 

behaviour of Gm and G M 2 decreases. Hence the approximation quality of G M 2 with G M 2 decreases, which results in 

a shape difference. 

To study the shift-orthogonal character of the ALO and LO pulses for various T, we have plotted the auto- 
correlations in Fig. [9] As can be seen, the samples r p T.°.M(mT) = p\ m as S m o, i.e. they vanish at almost each sam- 




(90) 
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Figure 6: ALO and LO pulses for Gaussian monocycle in time, L = 1 and M = 2K. 
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time in sec 



Figure 8: Centered ALO and LO pulse for T = 1.57" with L = 25, M = 2K = 40. 




Figure 9: Auto-correlation of pulses for L = 25, T = T p /K. 
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Figure 10: PSD of pulses for L = 25, T = T p /K. 



pie point except in the origin. The NESP performance for various values of T is shown in Fig. |T0j Approaching 
T = 7q cancels the FIR prefilter optimization of g, i.e. the spectrum becomes flat. Finally, in Fig.QT]the g am of our 
orthogonalization strategy can be seen. In both cases, the N-ary OOPPM design transmit at an uncoded bit rate of 



^.Sfio.Bw+i). (91) 

1 s 1 j\Jl o 

Fig.QT]shows the NESP value r\ over the transmit rate R/, for L = 25 and m = 2. Decreasing T = T p /K results in more 
overlap, which increases the number of symbols N in T s and hence Rt,, but only slightly decreases rj, see Fig. 1101111 

Summarizing, a triplication of the transmit rate from 0.18Gbit/s to 0.6Gbit/s is possible without loosing much 
signal power B/T s . Let us note the fact that the transmit bit data rate is an uncoded rate and is not an achievable rate. 
For ananalysis of achievable data rates by deriving the mutual information of the system see the work of Ramirez et. 
al in ifl^j and Giiney et. al in [62]. Obviously, the unshaped Gaussian monocycle then yields the highest transmit rate, 
since (|9TT i behaves logarithmically in the number of symbols, as seen in Fig.Qj] But this has practically zero SNR 
when respecting the FCC regulation and results in a high error rate (l87l).(l88l). On the other hand, a longer symbol 
duration, allows in © a higher energy S and hence a lower error rate in ( fTTb . 

Hence, the decreasing of the transmit rate due to the increased symbol duration used for FCC optimal filtering of 
the masked Gaussian monocycle can be more than compensated by the proposed OOPPM technique. 
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Figure 11: Orthogonalization results of the fixed FCC optimized pulse p for L = 25 in T s = 15GTo over various T resulting in various transmit 
rates R;,. 

7. Conclusion 

We have proposed a new pulse design method for UWB-IR which provides high spectral efficiency under FCC 
constraints and allowing a A^-ary OPPM transmission with finite transmission and receiving time by keeping almost 
orthogonality. In fact, the correlation parameters can be keep below the noise level by using small time-shifts T < T p . 
As a result, this provides much higher data rates as compared to BPSK or BPPM. Furthermore, our pulse design 
provides a N-ary orthogonal PSM transmission by getting a lower bit error rate at the price of a higher complexity. 

Simultaneous orthogonalization and spectral frequency shaping is a challenging and hard problem. We believe that 
for certain shifts being integer multiples of To, a numerical solver might be helpfully to directly solve the combined 
problem as discusses in Section 1531 

We highlight the broad application of the OOPPM design, not only being limited to UWB systems but rather 
applicable to a pulse shaped communication system under a local frequency constraint in general. 



Acknowledgements 

Thanks to Holger Boche for helpful discussions. This work was partly supported by the Deutsche Forschungsge- 
meinschaft (DFG) grants Bo 1734/13-1, WI 1044/25-1 (UKoLoS) and JU 2795/1-1. 



29 



Appendix A. 

Proof of Lemma\J] Let a, b € R with a + b = — 1. Then 5° and S h , defined by the spectral theorem, are also positive 
and self-adjoint on fi. Moreover for each / we have the following unique representation f - Yik c kPk with c e Z 2 (Z) 
due to the Riesz basis property. For / = p\ in (l72l i we get 

p, = J^WO* = Yj {SapuSbpk) -x Pk (A ' 1} 

k k 

since S"pi,S b pi l e W there exist unique sequences C/,dt s.t. S"pi = J^ a ci a p a ,S h pt : = Y^pdkpPp- Hence we get 



pi = ^ 2j Cla P a ' XI 



P* (A.2) 



= ^ c to^s(/'ff» (A.3) 

= 2 Z [C]to[G] ^ [D * ] ^ = 2 [CGD * ]ttW (A - 4) 

where c/„ and djg* are the coefficients of the biinfinite matrices C resp. D. Since for each / € Z we have Y^k $ikPk = Pi 
and{[CGD*] ft U<5 tt } el 2 , wegetO = £jt([CGD*] tt -6 lk ) Pk for each Ze Z. So by |Q Th.6.1.1(vuTl we can conclude 
that [CGD*]«. = 6 tk for all l,keZ and get 

CGD* = 1 o G = C^CD*) -1 o G 1 = D*C. (A.5) 

Obviously D and C are not an unique decomposition of S~ l , since a and b are not. If a — b — — |, we have D* = C* 
and hence G " = C = D. This establishes (l73l in an L 2 -sense. □ 

Appendix B. 

Proof of Theorem\3\ Let us first note that V(p) is a regular SI space since p is a stable generator. This has as con- 
sequence that frames are Riesz bases for V(p) ItSll Th. 2.2.7 (e)]. So any element / e V(p) = p *'c is uniquely 
determined by an i 1 sequence c. By the Riesz-Fischer Theorem this sequence c defines by its Fourier series a unique 
L 2 ([0, l])-function r = c. Hence, the Fourier transform of any / e V(p) is represented uniquely by r as / = rp, 
see also 11371 Th. 2. 10(d)]. On the other hand / is an orthonormal generator if and only if O/ = 1 a.e.. By using the 
periodicity of t we get 

*/ = X \P(- ~ k )\ 2 W - k )\ 2 (B-l) 

k 

= \T\ 2 Y J \P(--kf = \T\ 2 -® P = l. (B.2) 



Thus, we have |r| = 1/ y<t> p almost everywhere. Let us set f := 1/ y<t> ; , a.e. and a complex periodic phase function 
<p := e' a ^ : R — > {z e C | |z| = 1} with a : R — > [0, 27r] 1-periodic and measurable. Then any function t e L 2 ([Q, 1]) 
which satisfy ( 1B.21 > a.e. is given by t = f • <p a.e.. The L 2 -distance is then given by 
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tiP - f& = \\P - ft? = tiP - rPtili = \\P\\ 2 L 2 + \\rp\\ 2 L2 - f r\p\ 2 - f T\p\ 2 (B.3) 

JR JR 

= 2- f (t + t)|/3| 2 =2-2 f cos(ff) \p\ 2 ^ (B.4) 

JR JR 

> 2 -2 f |/3| 2 <V^ =2(1 - (B.5) 
Jk 

Since |/3| 2 is positive and <S> p is bounded and strictly positive a.e. the distance is minimized if and only if a(v) = 
a.e. in R, i.e. if we have equality in ( 1B.51 I. Hence <p = 1 a.e. and so f = t a.e., which corresponds hence to the unique 
orthonormal Lowdin generator / = p° with an L 2 -distance to p given in ( IB.5b . □ 

Remark. Note, that in fact the phase function has no influence on the power spectrum \p°,\ 2 - \(f>fp\ 2 - \fp\ 2 - \p°\ 2 . 
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